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We consider the effect of geometric frustration induced by the random distribution of loop lengths 
in the "fat" graphs of the dynamical triangulations model on coupled antiferromagnets. While the 
influence of such connectivity disorder is rather mild for ferromagnets in that an ordered phase 
persists and only the properties of the phase transition are substantially changed in some cases, 
any finite-temperature transition is wiped out due to frustration for some of the antiferromagnetic 
models. A wealth of different phenomena is observed: while for the annealed average of quantum 
gravity some graphs can adapt dynamically to allow the emergence of a Neel ordered phase, this is 
not possible for the quenched average, where a zero-temperature spin-glass phase appears instead. 
We relate the latter to the behaviour of conventional spin-glass models coupled to random graphs. 
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I. INTRODUCTION 

Models of random networks and surfaces have re- 
ceived extensive attention in statistical physics and field 
theory due to their wealth of applications in such di- 
verse fields as the modelling of the internet^, spin glass 
physics- and quantum gravity^. The diversity of appli- 
cations is reflected in a rather large variety of different 
graph models considered. Generic random graph and 
network models include the most general Erdos-Rcnyi 
model of n bonds distributed randomly between pairs 
selected from N vertices, as well as the Barabasi-Albert 
graphs constructed with the preferential attachment rule 
or the Watts-Strogatz method of interpolating between 
a regular lattice and a random graph via rewiring of 
links^. These constructions result in variable, but un- 
correlated vertex degrees with graph ensembles fully de- 
fined by the co-ordination number distribution P(q). On 
the other hand, fixed-degree or fc-regular random graphs 
have also been considered, which are equivalent to the 
non-planar or "thin" cf> 3 , 4 , ... Feynman diagrams of 
a zero-dimensional field theory^. While none of these 
networks feature a well-defined topology allowing for a 
local geometrical interpretation, a fattening of the thin 
graph propagators to ribbons yields orientable faces, en- 
abling an interpretation of the graphs as random sur- 
faces of fixed genus. The resulting dynamical triangu- 
lations model has been thoroughly studied with matrix 
model and combinatorial methods in the context of Eu- 
clidean quantum gravity^. The topological constraint in- 
duces spatial correlations in the degree distribution, and 
the ensemble of graphs is no longer fully determined by 
P{q)- Other ensembles of random graphs of well-defined 
topology have been studied, for instance the Poissonian 
Voronoi'-Delaunay tessellations resulting from a general- 
isation of the crystallographic Wigner-Seitz construction 
to randomly distributed generators^. 

A lot of effort has been invested in understanding the 
generic geometrical properties of these graph ensembles. 
In particular, for the network-type models dynamical 
properties have been regularly considered, such as the 



emergence of a giant component as edges are successively 
added to the graphsi. For the Dclaunay and dynamical 
tessellations with a well-defined topology, genuinely ge- 
ometrical attributes such as the fractal or Hausdorff di- 
mension and the correlation function of the local degrees 
have been of most interest^. As an alternative to direct 
investigations, these ensembles might be characterised by 
observing the cooperative behaviour of matter variables 
such as classical Ising spins si = ±1 with Hamiltonian 



J ■■ 



(1) 



or more general Potts or O(n) spins placed on the graph 
vertices, assuming ferromagnetic couplings Jy = Jo > 0. 
This exercise is also of interest from the inverted point-of- 
view of studying the effect of these various types of con- 
nectivity disorder on the spin models, being complemen- 
tary to the canonical case of weak disorder from random- 
ness of the couplings on regular lattices^. Graphs with- 
out topological constraint are inherently non-local and 
thus effectively infinite-dimensional. Consequently, one 
expects mean-field phase transitions for the ferromag- 
netic models. These are indeed found for the Erdos-Renyi 
and A;-regular (or thin) random gr aphs^^ i 8 1 9 1 1 , trans- 
forming them into convenient alternatives to treatments 
on the Bethe lattice or the complete graph, not encum- 
bered with boundary effects. More generally, degree dis- 
tributions with divergent moments can lead to interesting 
deviations from mean-field behaviour—. The effects on 
surface-like graphs are less homogeneous. For (uncor- 
related) bond disorder, the celebrated Harris criterion^ 
predicts a change of critical behaviour for cases with a 
positive specific-heat exponent a. An analogous criterion 
can be formulated for random graphs, taking into ac- 
count the spatial correlation of the degree distribution 1 ? , 
The predicted change in universality class for virtually 
all types of matter coupled to dynamical triangulations 
or fat graphs is in agreement with exact results for 
percolation^ and numerical investigations of the Ising 
and Potts models^. For Voronoi'-Delaunay tessellations, 
however, the predicted change of universality class for 
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FIG. 1: The link-flip move on two adjacent triangles of a 
dynamical triangulation (dashed lines) . The solid lines denote 
the corresponding dual, three-valent, "fat" <f> 3 graphs. 



models with a > such as the two-dimensional q = 3, 4 
Potts and three-dimensional Ising models is not observed 
numericall y 15 ' 16 . It is an open problem why, instead, 
these models behave according to their regular lattice 
critical exponents. 

In contrast to this weak-disorder case, the effect of con- 
nectivity disorder on models with antiferromagnetic in- 
teractions is much more profound. The existence of odd- 
length loops on many of the graphs discussed leads to 
severe geometric frustration with the possibility of al- 
together precluding the onset of a long-range ordered 
phase. Interestingly, however, this problem for random 
graphs has received very little attention to date^, such 
that the behaviour of, e.g., the Ising antiferromagnet for 
the various cases is unclear. The effect of tuning the 
amount of frustration on a regular lattice has been con- 
sidered in Ref. [Tt] for the ± J Ising model on the triangu- 
lar lattice, where it was found that a spin-glass phase ap- 
pears as the concentration p of antiferromagnetic bonds 
is increased from zero and disappears again as the sys- 
tem comes close to the perfect frustration of the pure 
antiferromagnet. Some use has been made of fc-regular 
or thin random graphs as an alternative realisation of 
the mean-field limit of spin glasses^, and in general one 
expects the perfect frustration of the effectively infinite- 
dimensional graphs of the non-topological type to lead 
to mean-field spin-glass behaviour for the antiferromag- 
net. The less extreme case of fixed-topology graphs is 
the subject of the present study. 

The rest of the paper is organised as follows. Sec- 
tion fll] discusses the general problem of frustration ex- 
erted by graphs of the dynamical triangulations type on 
antifcrromagncts and the annealed and quenched limits 
for the Ising antiferromagnet coupled to different graph 
types. The results of Monte Carlo simulations for the 
annealed case are discussed. In Sec. IIII1 wc investigate 
the quenched limit by means of exact ground-state com- 
putations in the framework of a defect-wall calculation, 
comparing Gaussian and bimodal spin glasses to the an- 
tiferromagnet. Finally, Sec. IIVI contains our conclusions. 



II. FRUSTRATION FROM FAT GRAPHS 

Dynamical triangulations originate in the discrete ap- 
proach to Euclidean quantum gravity via a path-integral 
quantisation of the gravitational interaction. The inte- 



gral over metrics is regularised by a sum over combina- 
torial manifolds or simplicial complexes 3 -. We focus here 
on the case of two dimensions, where this amounts to a 
sum over all possible gluings of equilateral triangles to 
closed, orientable surfaces of a given, for instance pla- 
nar or toroidal, topology. Graphs of this type are hence 
constructed by taking a number N2 of equilateral tri- 
angles and connecting them together along their edges 
at random, such that the resulting graph represents a 
closed and orientable surface of well-defined topology. In 
a computer simulation, such graphs can be implemented 
by the successive application of the "link-flip" move de- 
picted in Fig. [1] to the graph starting, e.g., from a regu- 
lar hexagonal lattice with periodic boundary conditions. 
The geometric duals of these triangulations naturally are 
graphs of fixed degree k = 3, known as planar 3 Feyn- 
man diagrams in field theory. The planarity constraint 
necessitates a "fattening" of the usual propagators to rib- 
bons, such that this type of diagrams is often refered to 
as "fat graphs"—. Instead of triangles, one can, of course, 
consider more general elementary polygons, such as the 
quadrangulations constructed from gluings of squares, 
whose duals are then consequently "fat" 4 graphs etc. 
The resulting combinatorial problem of counting these 
discrete surfaces can be solved exactly for various cases 
using a matrix-model formulation^ or alternative combi- 
natorial approache o 20 ' 21 . Geometrically, the most strik- 
ing result is that of an unusually large internal Hausdorff 
dimension dh = 4 for these topologically two-dimensional 
surfaces 3 -, resulting from a structure of "baby universes" 
connected to the main graph body with a minimal num- 
ber of links ("bottlenecks"). This fractal structure is ap- 
parent from the example graph embedding shown in Fig. 
[SJ while the graphs themselves are, by construction, de- 
fined in a purely intrinsic manner, without reference to 
an embedding into external space, an approximate em- 
bedding to some extent observing the constraint of equal 
edge lengths can be constructed to visualise the internal 
geometry^ 3 -. The distribution P(q) of co-ordination num- 
bers is found to fall off exponentiall y 2 ' 23 and hence does 
not exhibit the fat tails found in small-world networks. 
However, the structure of the graphs is not solely deter- 
mined by P(q), since the well-defined topology introduces 
long-range correlations in the co-ordination numbers, de- 
clining as r~ 2 at large separations 



A. Annealed average 

When coupling matter variables to the graphs, aver- 
ages have to be performed for the fluctuations of the 
graphs as well as those of the coupled spins with asso- 
ciated characteristic time scales T grap h and r ma ttcr- The 
annealed case of r grap h = T mat tor is the situation consid- 
ered in the context of matter coupled to quantum gravity. 
Since spin and disorder degrees-of-freedom fluctuate on 
the same time scales, the total free energy of the system 
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-Jo for all bonds, coupled to the following dynam- 



FIG. 2: Example of a dynamical triangulation of N2 = 5 000 
triangles embedded into Euclidean space. Note that the equi- 
laterality requirement for the triangles is only approximately 
fulfilled since an exact embedding into R 3 is not guaranteed 
to exist. 



is given by 



annealed 



(p) = -0- 1 ln(Z{P,g))c 



(2) 



where Z((3, Q) denotes the partition function for spins 
on a fixed graph Q and {-)g symbolises the average over 
the considered graph ensemble. Hence, the effect of ge- 
ometric randomness on the behaviour of spins is aug- 
mented by a back-reaction of matter variables on the un- 
derlying graphs. Formulating the combined problem as a 
multi-matrix model, several cases such as the (ferromag- 
netic) Ising2ii2£, g-state Potts-2£ZL and O(n) models 2 ^ 
(where the two latter classes, of course, include the sim- 
pler Ising model as the special cases q = 2 resp. n = 1) 
can be treated analytically (in the regime where they 
show continuous phase transitions on regular lattices, i.e., 
for q < 4 resp. —2 < n < 2). In these cases, the coupling 
to dynamical triangulations results in a shift of universal- 
ity class. More generally, in the framework of Liouville 
theory, the dressing of conformal weights A of critical 
matter coupled to two-dimensional quantum gravity is 
predicted to be2ML22 



A = 



VI -c + 24A- vT 



^/W~C- y/T 



(3) 



where c denotes the central charge^. The resulting crit- 
ical exponents agree with the exact results discussed 
above, and the predictions for a number of cases have 
been confirmed by numerical simulation studies, see, e.g., 

Refs. [mm. 

Coupling antiferromagnets to this type of random 
graphs leads to strong frustration^! due to the pres- 
ence of loops of odd length. Depending on the exact 
type of graphs considered, however, the annealed aver- 
age to some extent allows for the geometry to adapt to 
the antiferromagnetic interactions. In particular, con- 
sider an Ising antifcrromagnet according to Eq. {!]) with 



ical graphs: 

• Triangulations: Here, all elementary faces of the 
graph arc frustrated, and the frustration cannot be 
relieved by a dynamic response of the lattice. As 
has been shown by Wannier— , the Ising antifcr- 
romagnet on a regular triangular lattice remains 
paramagnetic down to zero temperature where, due 
to frustration, the system has a hnite residual en- 
tropy. It is clear that for any planar triangula- 
tion (whether regular or random) configurations of 
minimal energy have two satisfied and one bro- 
ken bond in each triangle. A large number of 
such states exists already for the triangular lattice, 
which is a member of the ensemble of dynamical 
triangulations considered. The freedom of chang- 
ing spin configurations without leaving the ground- 
state manifold is found to be local, offering no ener- 
getic reward and thus precluding long-range order. 
In fact, the residual, zero-temperature entropy of 
this model has been calculated^ to be Sq « 0.2613, 
to be compared to the value Sq ~ 0.3383 found for 
the triangular antifcrromagne1j2£ 

• Quadrangulations: Any planar quadrangulation 
is bipartite — it cannot contain any loops of odd 
length since all loops are composed of the elemen- 
tary faces of length four. This structure allows 
for a two-colouring of the lattice vertices in, say, 
black and white sites. Introducing extra (non- 
fluctuating) signs €i — ±1 at the vertices and per- 
forming a Mattis transformation as^£ 



Jij — Jij' 



(4) 



leaves the Hamiltonian ([T]) and thus the partition 
function invariant. Choosing = —1 for the 
white and ef = +1 for the black vertices maps 
the system identically to the Ising ferromagnet with 
Jij = +Jo- Hence, for the random lattice the re- 
sults of Refs. |24||25| apply, and the system under- 
goes a third-order phase transition to a Neel or- 
dered state. 

• <fi 3 or 4> 4 graphs: While neither of these graph 
classes is bipartite as a whole, it is clear that, e.g., 
with the hexagonal lattice (<fi 3 graphs) resp. the 
square lattice (</> 4 graphs) bipartite graphs occur in 
both ensembles. Hence the ground states of either 
system are perfectly Neel ordered configurations, 
and it is reasonable to expect long-range order to 
persist for some finite range of temperatures. 

To investigate this last case of a non-trivial interaction 
between geometric frustration and antiferromagnetism, 
we performed Monte Carlo simulations of the coupled 
system for graph sizes of N2 = 16 2 up to 128 2 vertices. 
As the temperature is lowered, we indeed find the emer- 
gence of antiferromagnetic Neel order with signatures of 
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FIG. 3: Magnetic susceptibility x = N 2 ({m 2 ) - (|m|> 2 ) of the 
Ising antiferromagnet coupled to annealed fat cj> 3 graphs with 
N 2 = 16 2 up to N 2 = 64 2 vertices. 



FIG. 4: Specific heat C v = N 2 p 2 {{e 2 } - (e) 2 ) of the Ising an- 
tiferromagnet coupled to annealed fat (f> 3 graphs from Monte 
Carlo simulations. 



a phase transition at rather low temperatures (3 C f=a 2.5, 
cf. the susceptibility data shown in Fig. [31 This tem- 
perature is to be compared with the phase transition at 
(3 C — | In w 0.773 found for the Ising ferromagnet on 
the "fat" <fi 3 graphs without self-energy and tadpole in- 
sertions considered here^. We find the specific heat to be 
completely independent of system size with a broad peak 
around (3 ss 1.4. cf. Fig. [4j and no crossing of the Binder 
parameter within the temperature range 0.1 < < 4.0 
considered. These findings, together with the apparent 
scaling of the susceptibility on the low-temperature side 
of the peak seen in Fig. [3j hint at the presence of a critical 
(low-temperature) phase with an associated Kostcrlitz- 
Thouless phase transition— This conjecture is further 
supported by the very slow convergence of effective criti- 
cal temperatures as the system size is increased, compati- 
ble with a logarithmic rather than a power-law approach. 
A series of simulations for the ferromagnet on the same 
lattices performed for comparison, on the other hand, 
yields (3 C = 0.783(7), d h v = 3.3(2) and y/d h v = 0.68(3), 
perfectly compatible with the exact result a 24 ' 25 . 

The mechanism of the antifcrromagnctic transition can 
be understood from an inspection of the distribution P(q) 
of loop lengths: as the temperature is lowered, the Ising 
antiferromagnet forces all odd- length loops out of the sys- 
tem, leading to surfaces completely composed out of faces 
of even length, i.e., squares, hexagons, octagons etc., cf. 
Fig. Thus, the back-reaction of the antifcrromagnctic 
matter on the graphs drives them into a non-frustrating 
phase of bipartite graphs compatible with Neel order. 
The resulting strong coupling between graphs and mat- 
ter throughout the whole low-temperature region could 
plausibly give rise to a critical phase with associated 
infinite-order phase transition as implied by the simu- 
lation results discussed above. This is in contrast to the 
Ising ferromagnet, where strong interaction between mat- 
ter and geometry is confined to the vicinity of the critical 
point. 



B. Quenched average 

The (graph-)qucnchcd limit T grap h ^ ^matter of aver- 
ages describes the case usually encountered in condensed 
matter physics, where impurities or lattice defects are 
fixed properties of the sample. Consequently, the average 
over disorder should be performed at the level of observ- 
able quantities, such as the free energy and its deriva- 
tives, leading to 

F q uc„chcd(/3) - -/3- 1 (lnZ(/?,g)) g , (5) 

This interchange of logarithm and disorder average as 
compared to the annealed limit of Eq. ([2]) often leads 
to dramatically different properties. Unfortunately, no 
exact prediction for ferromagnets in the spirit of ([3]) is 
available here. Since it is believed that ferromagnets with 
weak quenched disorder in two dimensions in general are 
related to conformal field theories with central charge 
c = 0, an approximation starting from the KPZ/DDK 
framework could be derived from Eq. ([3]) by setting c = 
there^ 2 -, i.e., 

A= -(V1 + 24A-1). (6) 

This form can also be more directly motivated by not- 
ing that the central charge of n replicas of matter vari- 
ables is additive, such that the limit n — * of the 
replica trick naturally leads to a central charge c = 
in (j3j 43 i 44 ] 45 . These predictions cannot be confronted 
with any exact solutions. Monte Carlo simulations for 
the Ising and Potts ferromagnets on quenched dynamical 
triangulation a 14 i 46 show a change in universality class, 
in agreement with the adapted relevance criterion for 
connectivity disorder discussed above^ 2 -. The critical ex- 
ponents, however, although changed from their regular- 
lattice values, in general seem not to be correctly pre- 
dicted by the form ([6]). 

For the antifcrromagncts, frustration from quenched, 
fat graphs is potentially stronger than in the annealed 
case: the frozen-in disorder of the graphs cannot adapt to 
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FIG. 5: Distribution P(q) of loop lengths of annealed, fat 
4> 3 random graphs coupled to an Ising antiferromagnet in the 
high-temperature phase at f3 = 0.1 and deep in the ordered 
phase at f3 — 3.2. Data for graphs with N2 = 32 2 vertices. 



alleviate the energy cost of odd-length loops. As before, 
the degree of frustration exerted on the Ising antiferro- 
magnet depends on the ensemble of graphs considered: 

• Triangulations: Following the argumentation 
presented above, the system is paramagnetic at 
all temperatures for each fixed random triangula- 
tion. Consequently, the quenched average merely 
describes the modulation of a completely disor- 
dered system. 

• Quadrangulations: The perfect bipartiteness of 
the graphs is not affected by the process of the dis- 
order average. Via the described Mattis transfor- 
mation, the system is identical to the Ising ferro- 
magnet on quenched quadrangulations. By univer- 
sality, the critical behaviour is expected to be that 
of the Ising ferromagnct on quenched fat 3 graphs 
studied in Ref . [hi 

• </> 3 or (p 4 graphs: This graph ensemble features 
a broad distribution of odd-length loops exerting 
strong frustration on the coupled antiferromagnet. 
The lack of a back-reaction of the Ising spins on 
the underlying graphs prevents the suppression of 
loops of odd length observed in the annealed limit. 
The bipartite graph configurations selected in the 
low-temperature phase of the latter only occur with 
vanishing weight in the quenched average. This 
appears to preclude the emergence of long-range 
order at finite temperatures. However, the frozen- 
in frustration might give rise to spin-glass order at 
zero temperature. 

The rest of this paper is devoted to an investigation 
of the possibility of such spin-glass order for the Ising 
antiferromagnet on quenched fat <fi 3 graphs. 



III. ZERO-TEMPERATURE PHASE ON 
QUENCHED FAT GRAPHS 

Due to the random distribution of faces of odd (frus- 
trated) and even (unfrustratcd) lengths, it is natural to 
suspect that the Ising antiferromagnet on 3 fat graphs 
behaves similarly to a ± J Edwards- Anderson spin glass 
on the same lattices. It is not a priori clear, how- 
ever, whether the long-range correlation of co-ordination 
numbers 1 - 2 , might cause any non-universal differences be- 
tween these cases. How can the appearance of spin glass 
be detected and its properties determined? Generalising 
Peierls' argument for the stability of an ordered phase, 
a droplet-scaling theory for the spin-glass phase can be 
formulated 4 ^. The role of the droplet surface (free) en- 
ergy is then taken on by the width J(L) of the distri- 
bution of random couplings for a real-space rcnormalisa- 
tion group decimation at length scale L. In the course of 
renormalisation, J(L) scales as J(L) ~ L e % defining the 
spin stiffness exponent 8 S . If the system scales to weak 
coupling, 8 S < 0, spin-glass order is unstable at finite 
temperature and the system is below its lower critical di- 
mension, with 9 S describing the properties of the critical 
point at temperature T = 0. On the other hand, 6 S > 
indicates stability of a spin-glass phase at finite temper- 
ature. Numerically, the domain-wall free energy can be 
determined from the energy difference between ground 
states of systems with different types of boundary condi- 
tions (BCs) chosen such as to induce a relative domain 
wall 4 !. 



A. Method and droplet length scale 

Following this programme, one requires the genera- 
tion of a number of statistically independent graph re- 
alisations for performing the quenched average. Subse- 
quently, for each realisation ground states of the Ising 
antiferromagnet should be calculated. An ergodic set of 
graph updates for dynamical triangulations (and the dual 
4> 3 graphs) is given by the so-called Pachner move s 23 ' 48 , 
which we employ to generate independent realisations of 
toroidal topology, starting out from a perfect hexago- 
nal lattice with periodic boundary conditions. Such a 
toroidal shape is needed to induce excited states by a 
change of boundary conditions as indicated above. The 
Monte Carlo equilibration is done for closed graphs, leav- 
ing the task of identifying appropriate loops along which 
to cut them open. Several possibilities come to mind: 

• In the original hexagonal lattice, the identification 
of boundaries is obvious. If, by construction, no 
link-flip updates are performed on the links crossing 
a boundary, the toroidal cuts are fixed in the equi- 
librated graphs. This simplicity, however, comes 
at the expense of strong boundary effects: due to 
the large fractal dimension 3 - dh = 4, each vertex is 
comparatively closer to one of the boundaries than 
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FIG. 6: On a planar graph, two interdigitating spanning trees 
on the original and dual lattices touch or cross all links. 

would be expected for a two-dimensional system. 

• Allowing flips of links crossing a boundary, one 
might keep track of the induced evolution of the 
boundary loops. This turns out to be rather in- 
volved, in particular due to the appearance of self- 
intersections of the meandering boundary lines. 
Consequently, we have not considered this ap- 
proach for the final ground-state computations. 

• Ideally, the graph generation should employ peri- 
odic boundary conditions preventing surface effects 
with the graphs being "cut open" once equilibrated. 
This amounts to the identification of topologically 
inequivalent elementary loops on the surface. It 
is not immediately obvious, however, how such an 
optimal homotopy basis could be computed for an 
intrinsic graph structure. 

As it turns out, a number of computational problems 
related to topologically motivated cuttings of surfaces are 
./VP-hard. In particular, identification of the so-called 
"cut graph" whose removal renders a surface topolog- 
ically trivial, i.e., the locus of all points with at least 
two shortest paths from a given base point, is a non- 
polynomial problem^. Rather surprisingly, however, a 
homotopy basis of minimum length in form of a system 
of loops joint in a common base point can be computed in 
polynomial time^. Such a system consists of 2g simple 
loops (where g is the genus of the surface), whose com- 
plement in the manifold is a topological disk. This de- 
composition proceeds with a variant of the tree-cotree de- 
composition proposed in Ref. Hi]: starting from a marked 
base point on a planar surface, the dual complement of 
a spanning tree is a spanning tree as well, cf. Fig. [5] For 
non-trivial topology, however, by Eulcr's formula both 
trees simply do not have enough edges to cover the whole 
graph, and it is easily seen that the edges not touched 
by either tree are those defining the topologically non- 
trivial loops with respect to the base point. A set of 
minimal loops can then be found by computing the tree 
T of shortest paths in the graph G and the maximum 
spanning tree T* in the complement (G\T)*, where the 
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FIG. 7: Scaling of the average length (L m i n ) of the domain 
wall in the ferromagnet with antiperiodic boundary condi- 
tions together with fits of the form (L m m) = AN^ dh ■ Note 
the considerably larger scaling corrections for the graphs with 
fixed boundaries. 



weight of each dual edge e* is chosen to be the length of 
the shortest loop in G containing the base point as well 
as the edge e. The desired basis is then given^ by the 
minimal loops defined by all edges e neither contained in 
T nor crossed by T*. 

Exact ground states of the Ising system on the re- 
sulting graphs arc computed in polynomial time via the 
mapping to a so-called minimum-weight perfect match- 
ing problem^. Due to a limitation of this approach only 
topologically planar graphs can be treated. Hence, after 
cutting the graph one of the resulting boundaries is left 
open, while cither periodic (P) or antiperiodic (AP) con- 
ditions are employed along the second boundary. The 
scaling of the energy difference, (\AE\p ij^>) ~ L^, 
then gives access to the spin-stiffness exponent 8 S . In 
contrast to the case of a regular lattice, where the rele- 
vant domain- wall length scale is simply equivalent to the 
length of the boundary, Low ~ L, the corresponding 
scale is not immediately obvious for the fractal graphs 
at hand. It can be found, however, from the ground- 
state problem of the ferromagnet: if AP boundary con- 
ditions are applied along one direction, the system re- 
sponds with a domain wall of minimum length L m i n , cor- 
responding to the length scale of the applied perturba- 
tion. The corresponding average length (L m i n ) is shown 
in Fig. for 5 000 graph replicas with N 2 < 256 2 . Here, 
(•) denotes an average over disorder. Fits of the form 

(L min ) = A L Nl /d ' h yield d' h = 4.0161(91) for the fixed 
boundaries resp. d' h = 4.0513(72) for the homotopy ba- 
sis, implying that in fact d' h — dh — 4, and the rele- 
vant domain- wall length scale is equivalent to the intrin- 
sic length L c ff ~ N^ dh of the fractal graphs, implying 

(|AJ5|p /A p) ~ N%' /dh . Note that the applied fits with- 
out correction terms allow an impressively precise deter- 
mination of dh, while it turned out to be very tedious to 
determine the Hausdorff dimension numerically by con- 
sidering the geometrical two-point functio n 23 i 53 . 
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spin glasses as well as the antiferromagnet on quenched, fat (j> 3 
random graphs. The lines show fits of the form (|A_E| P/ /ap) — 
E + S E N° s/dh to the data. 



formula ([3]) has been suggested for unfrustrated situa- 
tions, it is interesting to see what it predicts for the case 
at hand. Noting that 6 s /dh = —\jdh,v = A e — 1, a dress- 
ing according to ^ yields an invariant s /dh = for the 

±J model and a renormalised value 9 s /dh = —0.0886 
from the 9 s /d = —0.1422 found above for the Gaus- 
sian bonds. Our results of domain-wall energy compu- 
tations for (j> 3 random graphs of sizes N2 = 16 2 , 24 2 , 
. . ., 256 2 are collected in Fig. [5] together with fits of 
the form (\AE\ P/AP ) = E + S E N® s/dh . As for the 
honeycomb lattice, the Gaussian spin glass shows clear 
scaling (\AE\ p /ap) — > 0, with an estimated s /dh = 
—0.0684(25), different from the ordered lattice case, but 
also not quite compatible with the prediction of the 
quenched approximation (|6|). The ±J spin glass as 
well as the antiferromagnet show saturation behaviour 
implying 6 s /d h = with E = 0.599(23) (±J) and 
Eq = 0.472(44) (antiferromagnet), in agreement with the 
quenched KPZ prediction. In fact, also the finite-size ap- 
proach is very similar between these two cases, with es- 
sentially only the correction amplitude Se differing, cf. 
Fig. [5) This clearly supports the view that, on the pla- 
nar random graphs considered here, the antiferromagnet 
shows zero-temperature spin-glass behaviour in the uni- 
versality class of the ±J model, in contrast to regular 
lattices, where the antiferromagnet does not behave like 
a spin glass^. 



Domain-wall fractal dimension 



B. Spin-stiffness exponent 

To determine the spin-stiffness exponent 9 S , ground- 
state computations were performed for a series of graph 
sizes, comparing the ground-state energies for P and AP 
boundaries, using 5 000 disorder replica for each size. 
Since previous work has almost exclusively focused on the 
square lattice^ 4 -, as a benchmark and universality test we 
first considered the case of the regular honeycomb lattice, 
performing fits of the form {\AE\ P / AP ) = Eq + SeL 6b 
for L = 16, 24, . . ., 256. For a Gaussian distribution of 
the couplings Jij of ([I]), we arrive at Eq = 0.0026(261), 
clearly indicating that asymptotically Eq = as ex- 
pected. Fixing E = then yields 6 S = -0.2843(44) with 
quality-of-fit Q = 0.29. This is in excellent agreement 
with results for the square lattice^ 4 -, the negative value 
of 6 S indicating that spin-glass order is confined to zero 
temperature in this system. For symmetric, bimodal cou- 
plings = ± J, on the other hand, a clear saturation of 
defect energies is observed, resulting in Eq = 0.991(15), 
leading to a vanishing 6 S = of the effective spin-stiffness 
exponent. This is again in agreement with the square- 
lattice result, showing that the ± J spin glass is marginal 
in two dimensions (but it turns out that, in fact, T g = 
there^i, see also Refs. mil57|. 

Although the quenched approximation © to the KPZ 



It has been known for some time that domain walls 
in spin glasses are fractal curve o 54 i 58 . Recently, evidence 
has been presented that in two dimensions they arc in fact 
in agreement with the excursions known as "stochastic 
Loewner evolution" (SLE), which are closely related with 
unitary conformal field theorie a 59 i 60 . It is interesting to 
measure the fractal behaviour of domain walls on random 
graphs which are fractals themselves. Considering the 
link overlap qi = -^J2(i.j) s i s j s | s f f° r P/AP 
boundaries, it is clear that for Ising spins bonds crossed 
by a domain wall contribute —1 to qi, whereas all other 
bonds contribute +1 (Ni denotes the number of links of 
the graph). Hence one expects the scaling 

AT d 3 /d h 

where d s denotes the fractal dimension of the domain wall 
(note that = 3 N2 /2 from the Eulcr formula) . For the 
Gaussian spin glass on the honeycomb lattice, this ap- 
proach yields d s = 1.2725(33), in good agreement with 
the accepted value for the square lattic o 54 i 59 ' 61 . (Here, 
the smallest lattice sizes with L < 32 have been omit- 
ted from the fit to account for scaling corrections not 
explicitly taken into account.) For the bimodal distribu- 
tion, the comparison of P and AP ground states does not 
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FIG. 9: Scaling of zero-temperature domain-wall lengths for 
the Ising antiferromagnet as well as the ± J and Gaussian spin 
glass on fat 3 random graphs. The lines show fits of the form 
© to the data. 

define a unique domain wall due to the high degree of 
accidental degeneracy in the ground state, such that the 
above description randomly (but not necessarily without 
bias) captures one of these walls and a fit to the form 
([7]) yields dh = 1.283(11). Note that hence there is some 
uncertainty as to how to define and measure the fractal 
dimension properly in this case, and consideration of the 
backbone of ground states^, of domain walls of extremal 
length 61 or of the residual ground-state entropy 6 ^ lead to 
different estimates for the square-lattice system. Noting 
that d s /dh = 1-Adw, the quenched KPZ approximation 

([6]) predicts d s /dh — 0.4666 for the Gaussian spin glass 
coupled to the 4> 3 random graphs. The numerical results 
shown in Fig.[5]reveal clear scaling for this case with a re- 
sulting estimate d s /dh = 0.5042(24), again different from 
the regular lattice value d s /d = 0.63248(88) as well as 
from the quenched KPZ prediction. The corresponding 
estimates of d s /d h = 0.6330(10) and d s /d h = 0.6425(10) 
for the ± J and antiferromagnetic models due to degen- 
eracies again refer to some random choice of domain 
walls. More importantly, however, the scaling approaches 
for both cases arc again almost identical, cf. Fig. [5J sup- 
porting the view of a common universality class for both 



models. 



IV. CONCLUSIONS 

We have shown that the strongly frustrating influence 
of odd-length loops present in random graphs can trigger 
very interesting and characteristic effects in coupled an- 
tifcrromagncts, with the observed behaviour often being 
very different from that expected for fcrromagnets. Only 
for the special case of bipartite random lattices without 
odd-length loops such as dynamical quadrangulations, 
does a Mattis transformation identically relate the anti- 
ferromagnet to the ferromagnet. The annealed random- 
graph average considered in quantum gravity allows the 
lattice to dynamically alleviate frustration, leading, for 
instance, to the emergence of a Neel ordered phase for 
the (f> 3 and cf) 4 antiferromagnet, apparently accompanied 
by a Kosterlitz-Thouless phase transition (although, of 
course, it is not proven that this transition survives in the 
thermodynamic limit). Quenching the graphs, such dy- 
namical adaptation is precluded, moving the phase tran- 
sition to zero temperature, where instead a spin-glass 
phase appears. From an exact defect-wall calculation, it 
appears that this spin glass is in the same universality 
class as the ±J model on the same lattice. The wealth 
of these observed effects crucially depends on the local- 
ity (or finite dimensionality) of the considered surfaces, 
thus avoiding the possibly less interesting mean-field be- 
haviour induced by more generic random graphs. 
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